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Abstract 

We show that the thh'd-order optical response of disordered hnear J-aggregates can 
be calculated by considering only a limited number of transitions between (multi-) 
exciton states. We calculate the pump-probe absorption spectrum resulting from 
the truncated set of transitions and show that, apart from the blue wing of the 
induced absorption peak, it agrees well with the exact spectrum. 
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1 Introduction 



Since the discovery of the J-band of pseudoisocyanine (PIC) [1,2], the study 
of the collective optical properties of molecular aggregates has received much 
attention. The collective nature of the excitations gives rise to narrow ab- 
sorption lines (exchange narrowing) and ultra-fast spontaneous emission. The 
development of novel optical techniques, a continuously increasing theoretical 
understanding of the optical response of these systems, and the realization 
that this type of excitations play an important role in natural light harvesting 
systems have kept these materials in the spotlight. 

During the past 15 years a topic of particular interest has been the effect of 
multi-exciton states on nonlinear optical spectroscopies. Amongst the latter 
are the nonlinear absorption spectrum [3,4], photon echoes [5,6], pump-probe 
spectroscopy [6,7,8], and two-dimensional spectroscopy [9,10]. Prom a theoret- 
ical point of view, accounting for multi-exciton states requires the handling of 
large matrices, due to the extent of the associated Hilbert space. It is well- 
known, however, that for homogeneous linear aggregates, three states domi- 
nate the third-order response, namely the ground state, the lowest one-exciton 
state and the lowest two-exciton state. In practice, J-aggregates suffer from 
appreciable disorder, which leads to localization of the exciton states. It was 
shown in Ref. [7] that the three-state picture still holds to a good approx- 
imation, provided that one replaces the chain length by the typical exciton 
localization length. The existence of the so-called hidden structure of the Lif- 
shits tail of the density of states (DOS) justifies this approach [11,12]. 

In this paper, we put a firmer basis under the above idea, by systematically 
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analyzing the dominant ground-state to one-exciton transitions and one- to 
two-exciton transitions in disordered linear aggregates. By comparison to ex- 
act spectra, we show that, indeed, the pictme of three dominant states per 
localization segment holds. 



2 Model 



We model a single aggregate as a linear chain of coupled two-level monomers 
with parallel transition dipoles. We assume that the aggregate interacts with 
a disordered environment, resulting in random fluctuations in the molecular 
transition energies (diagonal disorder), and restrict ourselves to nearest 
neighbor excitation transfer interactions — J. The optical excitations are de- 
scribed by the Prenkel exciton Hamiltonian, 

N N-l 

H^J2^n\n){n\-Jj2{\n){n+l\+h.c.). (1) 

n=l n=l 

Here, \n) denotes the state with the nth site excited and all other sites in the 
ground state. The monomer excitation energies are modeled as uncorrelated 
Gaussian variables with zero mean and standard deviation a. For J-aggregates 
J > 0. Numerical diagonalization of the Hamilton yields the exciton energies 
Ejy {v — 1, . . . ,N) and exciton wavefunctions \u) — Y.n=i ^vn\n) of the one- 
exciton states, where ip„n is the nth component of the wavefunction The 
two exciton states are given by the Slater determinant of two different one- 
exciton states |z/i) and [z/a) [4]: |z/i,z/2) = Y.n<mi'P'^in<P'^2m - <Puim<Pu2n)\n,m) 
with |n, m) the state in which the sites n and m are excited and all other 
sites are in the ground state. The corresponding two-exciton energy is given 
by -Eivii/z = E^i + E^^. 
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3 Selecting dominant transitions 



The physical size of a molecular aggregate can amount to thousands of monomers, 
but the disorder localizes the exciton states [13]. For J-aggregates a small 
number of states at the bottom of the one-exciton band contain almost all the 
oscillator strength. The states are localized on segments with a typical exten- 
sion N*, called the locahzation length, which depends on the magnitude of the 
disorder [14]. The wavefunctions of these states overlap weakly and consist of 
(mainly) a single peak (they have no node within the localization segment), 
see Fig. 1. For the remainder of this paper we will refer to these states as s 
states, and S will denote the set of s states. 
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Fig. 1. 

Fig. 2. (a) The lowest 12 one-exciton states of a chain of length N = 500 for a 
particular disorder realization, (b) A subset of s states (black) and p states (gray) 
that mostly contribute to the one-to-two exciton transitions. 



From the complete set of wavefunctions we select the s states using the se- 
lection rule proposed in Ref. [12], | J2nVi^n\'^vn\ \ > C'o- For a disorder-free the 

lowest state contains 81% of the total oscillator strength between the ground 
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state and the one-exciton band [14]. Numerically we found that the s states 
selected by taking Co = 0.75 (for 0.05J < a < 0.3J) together contain 76% 
of the total oscillator strength: {J2ses A*s)/-^ ~ 0.76. Here /ig is the transition 
dipole moment from the ground state to the state In Fig. 3(a) we show 
the absorption spectrum calculated only with the s states and compare it to 
the exact spectrum. 

We observe that the s states give a good representation of this spectrum, 
except for its blue wing, where higher-energy exciton states, which contain 
one node within the localization segment, contribute as well [15,16]. These 
so-called p states may be identified as the second one-exciton state on a local- 
ization segment [11,12]. 

The p states play a crucial role in the third-order response. In order to analyze 
this, we have considered two-exciton states \s,p) given by the Slater determi- 
nant of a given s state (selected as described above) with all other one-exciton 
states V ^ S (the two-cxciton state consisting of two s type states localized 
on different segments do not contribute to the nonlinear response), and cal- 
culate the corresponding transition dipole moments Hsv,s- Ftom the whole set 
of iisi',s-i we select the largest one, denoted by iJ>sps,s, were the substrict s in 
Ps indicates its relation with the state |s). It turns out that the one-exciton 
state \ps) selected in this way is localized on the same segment as the state 
\s). Several of these doublets of s and p states are shown in Fig. 1(b). The 
partners of the lowest s states indeed look like p states, having a well-defined 
node within the localization segment. They form the hidden structure of the 
Lifshits tail of the DOS [11] we mentioned above. For higher lying s states, 
these partners (not shown) are more delocalized and often do not have a p-like 
shape. 
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The average ratio of the oscillator strength of the transitions |0) — > |s) and 
\s) \ps) turned out to be (^^%^^sl ^^t) ~ For the dominant ground-to- 
one and one-to-two exciton transition in a homogeneous chain, this ratio reads 
A*i2,i/a*i ~ 1-57. This comparison suggests that our selection of two-exciton 
states well captures the dominant one-to-two exciton transitions in disordered 
chains. We have also found that the energy separation between the states \s) 
and \ps) obeys {E^ - Ep^) ^ 1.35 x Sn^J/N*^, with N* = (//f/0.81). This 
resembles the level spacing of a homogeneous chain: E2 — Ei ^ 37i^J/N'^, 
confirming that the separation between the one-exciton bleaching peak and 
the one-to-two-exciton induced absorption peak may be used to extract the 
typical localization size from experiment [8]. 

To illustrate how the selected transitions, involving the {s,p) doublets, repro- 
duce the optical response of the aggregate, we have calculated the pump-probe 
spectrum at zero temperature using only these transitions [10], 

P{u) = ^ - 2f,tS{u;s -u;)+ f^XJiup^ - u;)), (2) 



and compared the result to the exact spectrum, see Fig. 3(b). Apart from a 
blue wing in the induced absorption part of the spectrum (the positive peak), 
the selected transitions reproduce the exact pump-probe spectrum very well. 
In particular the separation between the bleaching and induced absorption 
peaks is practically identical to the exact one. 
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Fig. 3. 

Fig. 4. (a) Absorption spectrum due to s states (circles) compared to the exact 
one (solid line), (b) Pump-probe spectrum due to the selected s and p type states 
(circles) compared to the exact one (solid line). The spectra were calculated for 
a = O.IJ (top) and a = 0.3 J (bottom). 



4 Conclusion 



We have shown that the third-order optical response of disordered linear J- 
aggregates is dominated by a very limited number of transitions. Considering 
only these transitions enormously reduces the computational effort necessary 
to simulate nonlinear experiments. We have used the procedure outlined above 
to calculate the optical bistable response of a thin film of J-aggregates, taking 
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into account the one-to-two exciton transitions [16] . 
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